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This paper addresses the optimal design and planning of the advanced hydrocarbon biofuel supply chain 
with the unit cost objective. Benefited from the drop-in properties of advanced hydrocarbon biofuels, 
the supply chain takes advantage of the existing petroleum infrastructure, which may lead to significant 
capital and transportation savings. A mixed-integer linear programming model is proposed to simul¬ 
taneously consider the supply chain design, integration strategy selection, and production planning. A 
robust optimization approach which tradeoffs the performance and conservatism is adopted to deal with 
the demand and supply uncertainty. Moreover, the unit cost objective makes the final products more 
cost-competitive. The resulting mixed-integer linear fractional programming model is solved by tailored 
optimization algorithm. County level cases in Illinois are analyzed and compared to show the advantage of 
the proposed optimization framework. The results show that the preconversion to petroleum-upgrading 
pathway is more economical when applying the unit cost objective. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Biofuels have been shown to be a promising fuel source of 
the future. They can be produced domestically from a wide vari¬ 
ety of biomass sources and reduce the dependence on fossil fuels 
(DOE, 2012). Moreover, greenhouse-gas emissions from biofuels 
are lower than those from their petroleum counterparts (Tong 
et al„ 2014a). Consequently, many countries have set national bio¬ 
fuels targets and provided incentives to accelerate the growth 
of bioenergy industry. In the U.S., the Renewable Fuels Standard 
(RFS), part of the Energy Independence and Security Act (EISA) 
of 2007, establishes an annual production target of 36 billion 
gallons of biofuels by 2022, of which 16 billion gallons should 
be advanced biofuels made from non-starch feedstocks to avoid 
adverse impacts on the food market (EISA, 2007). With the devel¬ 
opment of the third generation biofuel technologies, advanced 
biofuels can now be produced from cellulosic biomass such as 
crop residues, wood residues or dedicated energy crops (Yue 
et al„ 2014b). Moreover, advanced hydrocarbon biofuel products 
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(e.g. cellulosic-biomass-derived gasoline, diesel and aviation fuel) 
are functionally equivalent to the petroleum derivatives. Consid¬ 
ering all these promising properties, the expansion of the 
hydrocarbon biofuel industry is foreseeable in the next few 
decades, thus requiring the design and development of cost- 
effective biomass-to-biofuel supply chains (DOE, 2012). 

Many studies have been conducted on the design and plan¬ 
ning of biofuel supply chains from economic and environmental 
aspects (Akgul et al„ 2012; Bowling et al„ 2011; Dunnett et al„ 
2008; Elia et al„ 2011; Giarola et al., 2011; Sokhansanj et al., 2006; 
You et al., 2012; Yue and You, 2014; Yue et al., 2014a; Zamboni 
et al., 2009). However, the drop-in property of advanced hydrocar¬ 
bon biofuel is not well explored. The U.S. Department of Energy, 
DOE (2012) has pointed out the opportunities for the integra¬ 
tion of emerging hydrocarbon biofuel supply chains with existing 
petroleum production and distribution infrastructures. Although 
minimum retrofitting costs would be required in petroleum refin¬ 
ery for compatibility reasons, the integration indicates considerable 
capital savings on the construction of biofuel production facilities, 
which would help the biofuel products to be cost-competitive and 
bring in extra environmental benefits to the petroleum refineries 
(DOE, 2012). Researchers are investigating the possibility of con¬ 
verting biomass into biofuel in the traditional refinery. Huber and 
Corma (2007) summarized that catalytic cracking, hydrotreating, 
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and hydrocracking are the three main techniques for converting 
biomass into biofuel in existing petroleum refineries. DOE (2012) 
is investigating three possible insertion points into the petroleum 
refinery. They note that after converting the biomass into liquid 
bio-intermediates, they can be mixed with crude oil that feeds the 
Crude Distillation Units (CDU), or sent directly to upgrading units 
to produce gasoline and diesel. However, the studies on explicit 
economic evaluation of the integration possibility from the supply 
chain point of view are limited (Elia et al., 2012; Tong et al., 2014a, 
2014b). In this work, we propose and analyze a biofuel supply chain 
model integrating with existing petroleum refineries, combined 
with unit cost objective and the robust optimization framework. 

Most of the existing works consider the absolute economic per¬ 
formance, such as maximizing the total profit or minimizing the 
total cost. However, the economic objective along with per func¬ 
tional unit of the final product provides further opportunities for 
improving the economic performance, as the cost or profit would 
be reflected in the functioning outputs of the system: functional 
unit (Yue et al., 2013b). For instance, in the traditional supply chain 
models with cost minimization objective, the product demands are 
usually given with lower bounds. The optimal solution tends to 
produce and sell biofuels as less as possible to reduce the over¬ 
all cost. This is obviously not an economical solution, as the unit 
cost might be high. By using the unit objective, the optimal solu¬ 
tion can guarantee the lowest cost per unit, hence making the final 
products more cost-competitive. The unit objective can be defined 
as the total cost or profit divided by the total amount of functional 
unit. The problem can then be formulated as a mixed-integer linear 
fractional program (MILFP), which is a special type of non-convex 
mixed-integer nonlinear programs (MINLPs). Although the general 
MINLP solvers can be used, there might not be so efficient as some 
tailored solution algorithms. The parametric algorithm (Zhong and 
You, 2014) and reformulation-linearization approach (Yue et al., 
2013a) are two efficient tailored solution algorithms for MILFP 
problems as they take advantage of the efficient mixed-integer lin¬ 
ear programming (M1LP) methods to globally optimize the MILFP 
problems. Although unit objective has been chosen in many studies 
(You et al., 2009; Yue et al., 2013a, 2013b), the advantage of unit 
objective is rarely discussed. In this work, we present a detailed 
comparison between total cost minimization objective and unit 
cost minimization objective, and try to find out advantages of using 
the unit cost objective. 

The uncertainties in biofuel supply chain are critical and should 
be carefully considered. These uncertainties include seasonal and 
geographical fluctuation of biomass supply (Gebreslassie et al., 
2012; Tong et al., 2014a, 2014b), variability of biofuel demand 
due to unstable economic situations (Giarola et al., 2013; Kim 
et al., 2011; Marvin et al., 2012), fluctuating market price (Dal- 
Mas et al., 2011; Kim et al., 2011), and imprecise processing of 
data due to the process fluctuation and immature technologies 
(Tong et al„ 2014a). The inability to handle these uncertainties 
may lead to either an infeasible supply chain design or subopti- 
mal performance. The widely used approaches for optimization 
under uncertainty (Sahinidis, 2004) include stochastic program¬ 
ming (Birge and Louveaux, 1997; Tong et al., 2011), chance 
constraint programming (Charnes and Cooper, 1959; Yang et al., 
), and robust optimization ( i-Tal and Nemirovski, 2002; Li 
et al„ 2011; Li and Floudas, 2012). The goal of stochastic program¬ 
ming is to find the decision that is feasible for all the instances 
and maximizes the expectation of objective function over the ran¬ 
dom variables (Birge and Louveaux, 1997; McLean and Li, 2013). 
However it relies on the probability distributions of uncertain 
parameters. In the chance-constrained approach, uncertainties are 
represented through random variables with known probability dis¬ 
tribution and included in the constraints (Charnes and Cooper, 
1959). In robust optimization, only bounds of uncertain parameters 


are considered. Robust optimization aims to find the solution that 
is feasible under all the uncertainties (Ben-Tal et al„ 2009). Gen¬ 
erally speaking, a solution that is feasible for all the uncertainties 
usually does not leads to the best objective value, so the trade¬ 
off between robustness and performance is the main issue in the 
robust optimization. 

In this work, we address the optimal design and planning of 
advanced hydrocarbon biofuel supply chain under demand and 
supply uncertainties. Bertsimas and Sim’s robust counterpart opti¬ 
mization (Bertsimas and Sim, 2003; Li and Ierapetritou, 2008) is 
used to tradeoff the robustness and performance. Although uncer¬ 
tainties in biofuel supply chain optimization have been studied in 
many works, robust optimization is rarely used. A spatially explicit 
MILFP model with the unit cost objective is proposed to account 
for the drop-in properties of advance hydrocarbon biofuel supply 
chain, i.e. integrating with existing petroleum refinery infras¬ 
tructures. Moreover, two efficient algorithms, parametric method 
(Zhong and You, 2014) and reformulation-linearization (Yue et al., 
2013a), are adopted in this model, and their computational per¬ 
formance are compared with general purpose MINLP solvers. 
Additionally, the difference between total cost minimization and 
unit cost minimization are carefully analyzed and discussed based 
on the Illinois cases. The results show that the unit cost objective 
prefers the preconversion to petroleum-upgrading pathway to pro¬ 
duce biofuels. And the budget parameter f, which is defined as 
the maximal number of uncertain parameters that can reach their 
worst cases, provides a mechanism to tradeoff between conser¬ 
vatism and economic performance. 

The rest of this article is organized as follows. The background 
of biofuel supply chain and its major drop-in feature are high¬ 
lighted in the next section. This is followed by a formal problem 
statement. A brief introduction of the mathematical formulation of 
MILFP model is presented in Section 4. Detailed model formulation 
and nomenclature can be found in the supporting materials. The 
specified solution approaches are given in Section 5. A comprehen¬ 
sive comparison and analysis is presented in the case study section, 
and the concluding remarks are given at the end of the paper. 


2. Background 

The superstructure of the advanced hydrocarbon biofuel supply 
chain is illustrated in Fig. 1 . A typical advanced hydrocarbon biofuel 
supply chain consists of harvesting sites, biorefineries, preconver¬ 
sion facilities, upgrading facilities, and demand zones. The biomass 
is cultivated and harvested in harvesting sites. The harvested 
biomass feedstocks can either be sent to integrated biorefineries 
for direct production, or undergo a two-stage conversion process 
(Wright et al., 2008; You and Wang, 2011), namely preconver¬ 
sion and upgrading. Preconversion stage converts biomass into 
bio-intermediates (e.g. bio-oil and bio-slurry) that is economical 
and efficient for transportation, whereas upgrading stage upgrades 
the bio-intermediates into final products (e.g. gasoline, diesel) 
(Gebreslassie et al., 2013; Wang et al., 2013; Zhang et al., 2014). 

Both DOE (DOE, 2012; NABC, 2012) and UOP (Marker, 2005) 
noted that advanced hydrocarbon biofuel can take advantage of 
the existing petroleum refinery infrastructure, which means that 
it is possible to use the upgrading facilities in the existing refinery 
to upgrade bio-intermediates into final products. Moreover, most 
refineries in the United States do not operates at their full capac¬ 
ity, which provides a great opportunity for biofuel supply chains 
integration. DOE proposed three insertion points for integrating 
the biofuel supply chain with a petroleum refinery (DOE, 2012). 
In insertion point 1 (blue dashed lines), bio-slurry is first mixed 
with crude oil, then sent to CDU, and finally converted to gasoline, 
diesel and jet fuel through a series of upgrading units. In insertion 
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Fig. 1. Superstructure of advanced hydrocarbon I 
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point 2 (red dashed lines), bio-oil is directly sent to the upgrading 
units to produce biofuels. Note that in both insertion points 1 and 2, 
an additional pretreatment is needed to remove oxygen and other 
contaminants before the intermediates are sent to the operation 
units in refinery (Marker, 2005). For insertion point 3 (green dashed 
lines), biofuels can be blended with conventional fuels in petroleum 
refinery and then sent to customers using existing pipeline distri¬ 
bution network. 

3. Problem statement 

The problem we addressed in this work is stated as below. 

The superstructure of biofuel supply chain integrating with 
petroleum refineries are shown in Fig. 1. We are given a set of 
biomass harvesting sites, potential preconversion facility locations, 
possible upgrading facility locations, potential biorefinery loca¬ 
tions, existing petroleum refineries, and the demand zones. A set 
of biomass feedstocks (namely crop residues, wood residues, and 
energy crops) with their major properties, including moisture con¬ 
tent, harvesting cost, and the availability at each harvesting site 
are provided. The demand of fuels (gasoline, diesel, and jet fuel) 
at each demand zone is given with the upper and lower bounds. 
For each production facility (preconversion facility, upgrading facil¬ 
ity, and biorefinery), we are given a set of conversion technologies, 
capacity levels, and their corresponding conversion rates, operation 
costs, and capital costs. The data associated with insertion points 
in petroleum refinery is also given. These data include conversion 
rates, retrofitting capacity bounds, and capital costs. The available 
transportation modes, transportation capacities, and correspond¬ 
ing unit transportation costs and transportation distances are all 
known. We consider the annual case for our model, and the capital 
cost is annualized by a constant discount rate throughout the given 
project life time. 

The objective is to minimize the functional-unit-based eco¬ 
nomic performance in terms of unit cost, by determining two types 


of decisions: supply chain design decisions and supply chain plan¬ 
ning decisions. The supply chain design decisions include number, 
size, location, and technology selection of the preconversion facil¬ 
ities, upgrading facilities and the integrated biorefineries, and the 
insertion point determination for existing petroleum refineries. 
On the other hand, supply chain planning decisions are related 
to harvesting, production, and distribution. Harvesting decisions 
determine the biomass feedstock harvesting amount. Produc¬ 
tion planning includes feedstock consumption and product yield 
in preconversion facilities, upgrading facilities, biorefineries and 
petroleum refineries. The distribution decisions determine the 
transportation amount for each transportation link and mode. 

In our model, robust counterpart optimization approach is 
adopted to deal with uncertainties. The biomass availability at each 
harvesting site and the product demand at each demand zones 
are regarded as uncertain parameters. All the uncertainties are 
assumed to be symmetrically bounded, a e [a — a, a + 3], where a 
represents the uncertainty, a represents the nominal value, and & 
represents the variation amplitude. The budget parameter F, which 
defines the maximal number of uncertainties that can reach their 
worst case, is introduced to control the degree of the conservatism 
of the solution. 

4. Mathematical model formulation 

We develop an MILFP model addressing the optimal design and 
planning of advanced hydrocarbon biofuel supply chain integrat¬ 
ing with existing petroleum refineries with the unit cost objective. 
This model is a modification and simplification of the previous work 
by Tong et al. (2014a). The overall cost minimization objective is 
replaced by the unit cost minimization objective. The multiperiod 
model is replaced by the single period spatially explicit design 
model. Considering the length of the article, we give a brief intro¬ 
duction of this model in this section. A detailed description and the 
associated nomenclature are given in the supporting materials. 
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The objective is to minimize the functional-unit-based eco¬ 
nomic performance, which is defined as the annualized overall cost 
divided by the total biofuel in terms of Gasoline Gallon Equivalent 
(GGE) sold to the customers. The annualized cost involves the cap¬ 
ital cost for facility installation, operation cost for fuel production, 
and government incentives. 

m j n Qnnuaf _ Ccapital + Coperation ~ Cincentive , ^ , 

E P ,.,A ■ soId P,s/ ~ E P ,sA SOId P,sf 

where C annua i is the annualized overall cost, <p p is the coefficient for 
one unit of product p in volume to its energy content in terms of 
GGE, and sold p<s f is the amount of product p delivered to demand 
zone sf. Because both C annua i and sold ps f are variables, the non¬ 
convexity of the objective function leads to a nonconvex MILFP 
problem. 

The capital cost consists the cost from preconversion facilities 
tcapk k , upgrading facilities tcapn n , integrated biorefineries tcapli, 
and retrofit in petroleum refineries tcapsr sr , and is annualized by 
the discount rate IR. NY is the project life time. 

1R(1+/R) NY 


|yEcap/c fe + y^tcaprin + ^tcapli + y ^tcapsr sr j (2) 

The operation cost comes from biomass harvesting cost, opera¬ 
tion and maintenance (O & M) cost, variable production cost, and 
transportation cost. 

Coperation = Charvest + Q)M + C pr(H j Iicpon + Ctransport (3) 

The government incentives include the facility construction 
incentives and the volumetric incentives for biofuel production and 
usage (You et al., 2012). INCVOp is the unit volumetric incentive for 
biofuel p 

IR(\ +1R) NY 


and lerapetritou, 2008). There are several well-known robust coun¬ 
terpart optimization formulations, namely Soyster’s formulation 
(Soyster, 1973), Ben-Tal and Nemirovski’s formulation (Ben-Tal and 
Nemirovski, 1999) and Bertsimas and Sim’s formulation (Bertsimas 
and Sim, 2003), etc. However, the Soyster’s formulation is too 
conservative as it ensures the feasibility against all the poten¬ 
tial realizations, which sacrifices the optimal performance. The 
Ben-Tal and Nemirovski’s formulation provides a mechanism to 
allow a tradeoff between robustness and performance. However, 
the introduced nonlinearity complicates the computational perfor¬ 
mance. In our work, the Bertsimas and Sim’s approach is adopted, 
which not only avoids complicated nonlinear optimization, but also 
provides a way to consider the tradeoff between performance and 
conservatism (Bertsimas and Sim, 2003; Li and lerapetritou, 2008). 
In this robust counterpart optimization, a budget parameter r is 
introduced to control the degree of conservatism of the solution. It 
can take a value between 0 to the number of uncertain parameters 
in the constraints and not necessarily to be integer. It is unlikely that 
all the uncertain parameters get their worst values simultaneously, 
so the objective of this formulation is to use a user-defined value 
r to control up to parameters that are allowed to reach their 
worst case. The larger r is, the less likely constraint violation would 
be. 

The general optimization problem can be denoted as PI with 
the objective function (5) and constraint (6). 

(5) 

(6) 


PI : min cx 

s.t. y ai m x m < pi, V/ 


We assume that some of the parameter a lm are uncertain, and 
take values within the range of [a; m - a )m , a )m + a (m ], where a; m 
represents the nominal value and a ; m represents the variation 
amplitude. 

In Bertsimas and Sim’s robust formulation, a user-defined bud¬ 
get value ri is introduced to control up to [Fj parameters that are 
allowed to reach their worst case. Then the constraint (6) is first 
transformed into constraint (7) 


| y^inck k + y^incn„ + y^inc/, | + y'jNCVOp ■ sold p sf (4) 

The model satisfies the mass balance constraints, production 
constraints, and capacity constraints. The mass balance constraints 
define the mass conservation of materials at each node of the sup¬ 
ply chain. The production constraints describe the input-output 
relationship at the material processing facilities. Capacity con¬ 
straints limit the raw material purchases, production amounts, 
and the transportation flows. In our model, the drop-in property 
in advanced hydrocarbon biofuel supply chain (integrating with 
existing petroleum refineries) can be appropriately modeled by 
using the binary variables isrl sr and isr2 sr , which denote whether 
petroleum refinery sr is retrofitted by insertion point 1 or insertion 
point 2. The readers can refer to the supporting information for the 
detailed description of the proposed model formulation. 

5. Solution approach 

5.1. Robust optimization 

Robust optimization usually outperforms other methods to 
handle uncertainties due to its independence on the probability 
distribution of uncertain parameters. The objective of the robust 
counterpart optimization is to choose a solution that is able to 
cope best with the various realizations of the uncertainty data (Li 


y]a /m x m + max jy^a /m |x m | + (/]- [T,J) • a /t| |x t ,11 < P/, VI 

(7) 

where S; represents the subset that contains I TjJ uncertain param¬ 
eters, and t; is an index to describe an additional uncertain 
parameter if r t is not an integer. In order to remove the abso¬ 
lute value |x m ], an additional nonnegative variable u m = \x m \ is 
introduced. Then the constraint (7) can be reformulated as 

J2 a lm x m + max | ya Im u m + (r, - [T,J) • a ltf u t , j < p h VI 


Note that the maximization problem in (8) is equivalent to prob¬ 
lem (10)—(12), because the optimal solutionz ( * m must consist of [/j J 
variables at 1 and one variable at /j - [/jj. 


max y]a, m u m z, m 
s.t. y Zlm <r h VI 


(10) 

(11) 
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0 < Z; m < 1, VI, m (12) 

The dual problem of (10)-(12) is described as follows 
min nzi+^im (13) 

meM, 

s.t. z, + q lm >a lm u m , VI, m (14) 


Similarly, we assume the biomass availability BA bi takes value 
in the range of [BA bi - BA b h BA b t + BA b Then after introducing 
the budget parameter r b i with the range of [0,1 ], the original con¬ 
straint (30) can be reformulated into (31) and (32). 

harv b i < BA hh Vi, b (30) 

harv b i - BA b i + r M • z M + q bJ < 0, Vb,i (31) 


q, m >0, VI, m (15) 

z, > 0, VI (16) 

where q\ m corresponds to the dual variable of the inequality Z\ m < 1, 
Z\ corresponds to the dual variable of inequality E m jZ im < ri- 
The robust formulation can then be transformed to the follow¬ 
ing counterpart P2 by substituting the inner optimization problem 
with the equivalent optimization problem. 


P2 : min cx (17) 

s.t. ^a^Xm + r^+^qta^p,, VI (18) 

m meMi 

Zi + qim>Oi m u m , VI, m (19) 

-Um<Xm<U m , Vm (20) 

qim> 0, VI, m (21) 

z, > 0, VI, m (22) 


For this robust counterpart, the probability bounds of con¬ 
straints violation can be calculated by the number of uncertain 
parameters n, and the budget parameter r b For more details about 
this probability bounds calculation, please refer to the work by Li 
and Ierapetritou (2008) and Bertsimas and Sim (2003). 

In our work, fuel demands and biomass availabilities are 
considered uncertain. The demand uncertainties are denoted 
as DFMp^ 18 , which take value in the symmetric region 

of DEM“% lb e [DEM u p p ' LB - DEM u p % LB , DEM^f + DEM^j . The 
original demand constraints are listed by Eqs. (23) and (24). 

S0ld p sf > DEM LB sf Vp, sf (23) 


Zb,i + qb,i>BA b 'i, Vb,i (32) 

5.2. Mixed-integer linear fractional programming 

Note that in our model, only the objective function has the 
linear fractional form, whereas other constraints are all linear. It 
can be regarded as an MILFP problem, which is a class of non- 
convex MINLP problems. The general-purpose solvers for MINLP 
problems are DICOPT, SBB, and BARON, etc. DICOPT is based 
on the Outer-Approximation algorithm, but assumes the model 
to be convex (Duran and Grossmann, 1986). Applying DICOPT 
algorithm in the non-convex optimization problems may lead 
to suboptimal solutions. SBB implements the simple branch and 
bound method (Bussieck and Drud, 2001 ). Also, it only guarantees 
global optimality for convex MINLP problems. BARON is a global 
optimization algorithm and based on the branch-and-reduce algo¬ 
rithm (Sahinidis, 1996). The computational performances of these 
general-purpose solvers are not satisfactory, because they do not 
take advantage of the special structure of MILFP problems. In this 
work, we applied two tailored solution approaches for MILFP prob¬ 
lems, namely parametric approach (Zhong and You, 2014) and 
reformulation-linearization approach (Yue et al., 2013a), and com¬ 
pare the performance with the MINLP solvers. 

min ^o + Ei^ + E/ 2 ^ 

B0 + Ei B1 iXi + Ej B2 jXi (33) 

s.t. C0 fc + ^Cl ffc Xi + = 0, Vk 


The MILFP model in this paper can be generalized by P3. First 
we introduce the parametric approach. The main idea is to itera¬ 
tively solve a series of equivalent MILP subproblems. These MILP 
problems are given by problem P4. 


S0ld p sf < DEM p ’’ f Vp, sf (24) 

The uncertainties appears in the right hand side of the con¬ 
straints. To apply the robust formulation to the demand constraint 
(23), an additional variable with fixed value of 1 is added and the 
constraints is formulated as follows: 


min ^AO + ^Al*+ ^ 

s.t. COfc + y^ClftXj + y^C20j = 0, Vk 


E B i' Xi+ E B2 ^) 

(34) 


-SOldprf + DEM^ B sf • 1 < 0, Vp, Sf (25) 

Thus, this constraint has only one uncertain coefficient, and the 
budget parameter r 18 ^ takes a value in the range of [0,1]. The fol¬ 
lowing constraints (26) and (27) are then incorporated into the 


robust optimization. 

-SOldp^f + DEM p B sf + z“ s/ • r PB f + q L p B sf < 0, Vp, sf (26) 

z p B sf + <?“/ ^ > v P>s/ (27) 

Similarly, for the demand upper bound constraints, we can get 
the following constraints 

soldp sf - DEM p p sf + r p 8 f ■ z p 8 f + q p p f < 0, Vp,sf + (28) 

z p!sf + 9p B / — DEM p p f , Vp, sf (29) 


The subproblem P4 has exactly the same constraints as the original 
MILFP problem P3, but with a linear objective function, instead of 
a nonlinear one. Following the idea of Newton’s method, the main 
steps for the algorithm can be summarized as following. 

Step 1. By setting r= 1, initialize q r to 0. 

Step 2. Solve PI and denote the optimal solution as x r and y r . 

Step 3. If the optimal value in PI is less than the optimality toler¬ 
ance 8, stop and output (x r ,y r ) as the optimal solution; otherwise, 

A0+E Al i*. +E A2,yT 

letq r+ i = B0+ sp B1 .^r + Y^ m y • g° to step 2 and replace q r with q r+1 . 

Another approach is the reformulation-linearization method 
(Yue et al., 2013a). The main idea is to reformulate the MILFP 
problem into an exactly equivalent MILP problem by introducing 
additional continuous variables and constraints. 
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Firstly, nonnegative variables u and Xj, such that 
u = —-—\—'—-——-—, and Xj = —— . . —-——-— = u • Xj are 
B0 +E, B1 < x i + Ej B2 J y j B0 +E, B1 < x i + Ej B2 j y j 

introduced. Then the objective function in P3 can be reformulated 
as min AO - u + ; A 1 jX, + Ei^ 2 i Y i< where the additional variables 

Yj = u yj is introduced to remove nonlinearity. By multiplying both 
sides of the constraints in P3 by u, the resulting constraints can be 
reformulated as CO k u + X; + JjC2 ji( ’ Y i = °> vk - Moreover, 

the only nonlinear term Yj = uyj can be linearized by using the 
Glover’s linearization scheme (Glover, 1975), where the Big-M 
is introduced. Therefore, the MILFP problem P3 can be exactly 
reformulated to an M1LP problem P5. 

P5 : min AO - u + ^AljX; + ^A2jYj 
i j 

s.t. CO k u + ^Cl ik • Xj + ^C2 jfc • Yj = 0, Vk 

B0 u + ^Blj • X, + Y^ B2 i ■ Y i = 1 (35) 

Yj < u, Vj 
Yj < M-yj, Vj 
Y] >u~M.[l-yj), Vj 

There is exactly one to one mapping between the original vari¬ 
ables Xj,yj and the introduced variables Xj, Yj. Continuous variable 
x,- do not exist in the new formulation. Problem P5 only needs to 
be solved once, but it is more computationally expensive as some 
extra variables and constraints are introduced. These properties are 
reflected and carefully analyzed in the case study section. 


5.3. Integrated optimization framework for robust optimization 
and MILFP 

In this section, the MILFP and robust counterpart optimiza¬ 
tion are integrated in our optimization framework. The solution 
strategy of this integrated model is depicted as follows. When 
considering the unit cost objective, we transform the original sup¬ 
ply chain model into MILFP problem by changing the total cost 
minimization objective into the unit cost minimization one (1). 
Then robust counterpart programming can be applied into the 
MILFP model to deal with the demand and supply uncertain¬ 
ties. Note that when applying the robust optimization approach, 
only the constraints with uncertain parameters are modified and 
some additional constraints are introduced. Therefore, the result¬ 
ing robust counterpart formulation is still MILFP model, and can 
be solved by either general-purpose MINLP solvers or the tailored 
approaches presented in this paper. These steps are summarized in 
Fig. 2. 





6. Case study 

6.1. Input data 

To illustrate the application of the proposed framework, we 
solved a series of county level cases for the state of Illinois. Except 
the data related to time periods, all the other data are the same 
as the one in the previous work (Tong et al., 2014a), which are 
obtained from technical reports (DOE, 2011, 2012) and existing 
literatures (You and Wang, 2011 ). 

The state of Illinois is comprised of 102 counties, each of which 
is recognized as a node in our model that represents a harvesting 
site, a potential preconversion facility, a possible upgrading center, 
a potential integrated biorefinery, or a demand zone. Three major 
types of biomass are considered in our cases: crop residues (e.g. 
corn stover), energy crops (e.g. switchgrass), and wood residues 
(e.g. forest residues). We consider three types of liquid fuel prod¬ 
ucts, namely gasoline, diesel and jet fuel. Their demand are assumed 
to be proportional to the population in each county. In our cases, 
we consider two integrated conversion pathways for the integrated 
biorefinery (gasification + FT synthesis and pyrolysis followed by 
hydroprocessing) (Swanson et al., 2010; Wright et al., 2010), two 
preconversion technologies (rotating cone reactor pyrolysis and 
fluidized bed reactor pyrolysis) (Magalhaes et al., 2009; Uslu et al„ 
2008), and three upgrading technologies (bio oil to fuel using gasi¬ 
fication, bioslurry to fuel using gasification, and bio oil to fuel using 
hydroprocessing) (Gebreslassie et al., 2013; Magalhaes et al„ 2009; 
Uslu et al., 2008 ). Three capacity levels are considered for all the bio¬ 
fuel production facilities (preconversion facility, upgrading facility, 
and biorefinery). 

There are four existing petroleum refineries located in Illinois, 
namely the Lemont, Joliet, Robinson, and Wood River refineries. 
Their locations and capacity levels are known. We assume that 
each petroleum refinery operates at 95% of its capacity level for 
processing crude oil, while the remaining 5% capacity could be used 
for upgrading bio-intermediate (Tong et al., 2014a). In addition, 
three insertion points are analyzed in our case. 

In order to reduce the complexity of the supply chain optimiza¬ 
tion problem, we use the same heuristic model reduction rules as 
in our previous work (Tong et al., 2014a) to reduce candidate sites. 
After reduction, the supply chain network contains 40 harvesting 
sites, 19 possible preconversion facilities, 20 potential biorefineries, 
20 possible upgrading facilities, and 102 demand zones. 

All the cases were performed on a Dell desktop with an Intel 
Core i7-930 2.80 GHz and 6 GB RAM. All the models are coded in 
the GAMS 24.1.3 (Rosenthal, 2012). The M1LP model is solved with 
the solver CPLEX 12.5, and the MINLP model is solved by DICOPT, 
SBB, BARON 12.7 solvers. The optimal gap in all cases is set to 0. A 
time limit of 10 h is applied to all the models. 

6.2. Mixed-integer fractional programming 

In this subsection, we present a comprehensive comparison 
between the supply chain optimization problems with total cost 
minimization objective and those with unit cost minimization 
objective. 

We first define two models considered in this subsection, the 
deterministic model and the MILFP model. In deterministic model, 
all the parameters are deterministic. Moreover, the objective is to 
minimize the annualized cost. 

min C annual = C ca pj tal + C operation -C incent j Ve (36) 

The deterministic model consists of the objective function (36) 
and the constraints (SI )-(S52), (S54)-(S80) in the supporting infor¬ 
mation. In contract to the deterministic model, we use the unit cost 
as the objective function in the MILFP model. The MILFP model 
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consists of the objective function (1) and the constraints same as 
those in deterministic model. This model can be solved by either 
MINLP solvers, such as D1COPT, SBB, BARON, or the parametric and 
reformulation-linearization approaches. 

First, we compare the economic and computational perfor¬ 
mance between the deterministic model and MILFP model. The 
computational result is shown in Table 1 . The total annualized cost 



is $2229 million in the deterministic model, which is 8% lower 
than the one in the MILFP models ($2421 million by parametric 
approach). However, the unit cost is $3,180 per GGE, which is much 
higher than the one in the MILFP models. This shows that mini¬ 
mizing the total annualized cost may increase unit cost to some 
extent. We also compare the computational results between the 
tailored algorithms with the general-purpose MINLP solvers. The 


Table 1 

Computational performance in deterministic model and MILFP model. 


Model and algorithm 


Total cost (MM $) Unit cost ($/GGE) a 


Deterministic 2229 

Parametric 2421 

Reformulation 11 2424 

MILFP model DICOPT 2421 

SBB 2573 

BARON 12.7 3461 


3.180 

2.979 

2.981 (2.873) 
2.979 

3.180(2.654) 
4.223 (2.699) 


a With the lower bound in the parentheses. 
b Reformulation means reformulation-linearization method. 
c Out of memory error. 

d Computational time limit of 10 h is applied to all 


Time(s) Gap No. of binary 


169 0% 422 

889 0% 422 

17,052 c 3.65% 560 

3577 0% 422 

36,000 d 16.53% 422 

36,000 d 36.09% 422 


No. of constraints 


32.731 
32,734 
40,687 

32.732 
32,732 
32,732 


I the algorithms. 
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Fig. 5. Cost breakdown in four models. 


parametric method gets the optimal solution of $2.979/GGE within 
889 s, which is demonstrated to be the most efficient approach. 
The reformulation-linearization approach terminates at the unit 
cost of $2.981/GGE at around 4.74 h when it encounters the out 
of memory error. Although it stopped before time limit and it has 
a large optimality gap (3.65%), the unit cost is quite close to the 
optimal result ($2,981 /GGE vs. S2.979/GGE). DICOPT solver returns 
the same optimal solution as parametric approach within 3577 s, 
although it can only guarantee the optimal solution for convex 
MINLP problems. By contrast, the computational performances of 
SBB and BARON 12.7 are not satisfactory. They stop at 10 h time 
limit with the optimality gap of 16.53% and 36.09%, respectively. 
Their optimal values are also far away from the global optimal 
value. Note that in MILFP model, the ones using reformulation- 
linearization approach, SBB solver, and BARON 12.7 solver do not 
return the zero gap at the 10 h time limit. We tested MILFP model 
using BARON and SBB solvers with lOOh time limit, and try to 
find more information. Unfortunately, the optimal solutions in both 
cases do not increase in the last 90 h. The lower bound in BARON 
increases from S2.692/GGE to S2.728/GGE, with the optimality gap 
decreases from 36.25% to 35.40%. The lower bound in SBB increases 
from S2.654/GGE to S2.668/GGE, with the gap decreasing from 
16.54% to 16.10%. This shows that increasing computational time 
does not improve the solution too much, and 10 h time limit is 
acceptable for these cases. The detailed comparison between 10 h 
and 100 h computational time limit is provided in the supporting 
information. 

The optimal design for the deterministic model is shown in 
Fig. 3, whereas the optimal design for MILFP model obtained from 
parametric approach is shown in Fig. 4. We can see several signif¬ 
icant differences between these two models. In the deterministic 
model, 10 preconversion facilities are built with the capacities ran¬ 
ging from 330 to 894kton per year. The preconversion facilities 
in Iroquois and Livingston Counties use the fluidized bed reactor 
pyrolysis technology, which produces the bio-slurry, while oth¬ 
ers use the rotating cone reactor pyrolysis technology to produce 


bio-oil. Eight biorefineries are built with the capacities ranging from 
17 to 91 MM GGE per year, all of which use fast pyrolysis followed 
by hydroprocessing. In the MILFP model, 12 preconversion facilities 
are constructed with the capacity ranging from 235 to 1331 kton 
per year. Four preconversion facilities using fluidized bed reactor 
are built in Christian, Pike, Vermilion, and Wayne Counties. Others 
use the rotating cone reactor technology. Six biorefineries are built 
with capacities ranging from 16 to 30 MM GGE per year, all of which 
use fast pyrolysis followed by hydroprocessing. Additionally, one 
upgrading facility with technology of hydroprocessing converting 
bio-oil to liquids is built in Kane County. Apparently, the unit cost 
in the MILFP model is lower than the one in deterministic model, 
whereas the overall cost is higher. And we have more preconversion 
facilities and fewer biorefineries in the MILFP model than those in 
deterministic model. In the MILFP model, more fuels are produced, 
especially in petroleum refinery. This is because, in the determi¬ 
nistic model with overall cost minimization objective, the decision 
makers tends to meet the minimal demand requirement in order to 
minimize the production cost, transportation cost, harvesting cost, 
etc. However, in the MILFP model, we tend to minimize the unit cost 
per GGE sold to the customers. It is more related to profitability to 
some extent. And it is a tradeoff between the overall cost and total 
amount of product sold. From the unit cost point of view, the pre¬ 
conversion to petroleum upgrading pathway is more economical, 
which leads to more production in refinery. 

The cost breakdown information in these two models is shown 
in Fig. 5 (a for the deterministic model and b for the MILFP model). 
We can see in the MILFP model, the proportion of capital cost, 
decreases from 16% to 13% compared with deterministic model. In 
the meantime, other costs, such as transportation cost, harvesting 
cost, and O & M cost increase slightly. It can be concluded that we 
can make unit cost lower by changing the supply chain structure at 
the expense of increasing transportation cost and production cost. 

As can be seen from Fig. 6 and Table 2, all these four refiner¬ 
ies are retrofitted. However, in the deterministic model, only the 
Joliet refinery used all the remaining capacity in refinery (5%). The 
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rofit information 


in four models. 


Refinery 


MILFP 


l a 2 % b 1 2 % 


2 % 1 2 % 



Wood River 

Total 129 


95 

53 

70 

219 


80% 87 

100% 154 

47% 108 

42% 21 147 

67% 129 388 


100% 76 

100% 49 106 

96% 51 

100% 75 

99% 124 232 


a 1 and 2 stands for the retrofitted capacity (MM GGE per year) in refinery by insertion 1 and 2, respectively. 
b % stands for the percentage of unused capacity retrofitted in refineries. 


26 


46% 70 

45% 29 

68% 124 


100% 

100% 

100% 

100% 

100% 





Fig. 6. Refienry retrofit information in four models. 


overall unused refinery capacity utilization is 67%. This is the best 
configuration in the cost minimization problem. Upgrading bio¬ 
fuel in petroleum refinery is a good alternative, but not always the 
best way to produce biofuel. However, in the MILFP model, almost 
all the petroleum refineries use their unused capacity to upgrade 
biofuels. Moreover, one upgrading facility is built instead of using 
biorefineries. It is a more economical way in terms of unit cost. 

Fig. 7 shows the important market share in these two models. 
Market share in this paper is defined as the percentage of biofu¬ 
els produced in different facilities (biorefinery, upgrading facility, 
and petroleum refinery). The difference of market share between 
the deterministic model and the MILFP model is significant. First 
of all, the products upgraded in petroleum refineries in the MILFP 
model is greater than those in the deterministic model (538 MM 
GGE per year compared to 364 MM GGE per year). Moreover, in 
the deterministic model, the petroleum refineries produce 52.5% 
of total biofuels, whereas biorefineries produces the other 47.5%. 



While in the MILFP model, petroleum refineries produce 66.1 % bio¬ 
fuels, biorefineries 16.8%, and upgrading facilities 17.1%. This again 
shows that the pathway of preconversion to petroleum upgrad¬ 
ing is more economic from the unit cost point of view. Moreover, 
in the MILFP model, all the unused capacity in petroleum refiner¬ 
ies are utilized for upgrading biofuels, and an additional upgrading 
facility is built to upgrade biofuels instead of producing them in the 
integrated biorefineries. 

6.3. Robust optimization 

In this subsection, another two models are introduced, the 
robust model and the robust-MILFP model. Robust model consists 
of the objective function (36), the robust counterpart (26)-(29), 
(31 )-(32), and other constraints (S2)-(S51), (S54)-(S80) in the sup¬ 
porting information. 

We first compare the results for the robust model with differ¬ 
ent budget value r in order to choose a good budget value for 
robust model in comparison. We set the demand uncertainty with 
5% variability, biomass availability with 10% variability. Then we 
solve this problem by robust optimization with 5 different budget 
parameters r ranging from 0 to 1. As noted previously, the larger 
r is, the more conservative the result would be. Therefore, both 
the annualized cost and unit cost increase as the r increases. The 
computational results listed in Table 3 also support this conclusion. 
In Fig. 8, the increase of budget parameter, which means less vio¬ 
lation of constraints, contributes to the increase of transportation 
cost, harvesting cost, and production cost. However, the O 8; M cost 
and capital cost almost remain constant. 

We choose a budget value of 0.75 in our robust model for 
later comparison with the deterministic model, which is a tradeoff 
between the conservatism and economic performance. The com¬ 
putational results are shown in Table 4. In this robust model, the 
annualized cost of $2354 million and the unit cost of $3,244 per 
GGE are a little higher than those in the deterministic model due to 
the conservatism of the robust formulation. The computation time 
is 109 CPUs. It shows that Bertsimas and $im’s formulation do not 
increase the computational time. 


1000 



■ biorefinery ■ petroleur 


■ upgrading 


■T=0 ^=0.25 


■r=0.75 «r=i 


Fig. 7. Market share of production facilities in four models. 


Fig. 8. Cost breakdown in robust model with different budget value V. 
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Table 3 

Computational result in robust model with different budget value r. 


r Transportation (MM $) 0&M(MM$) Production (MM $) Harvest (MM $) Capital (MM $) Annual (MM $) Unit cost (S/GGE) Time(s) 


0 402.9 

0.25 431.8 

0.5 441.4 

0.75 455.5 

1 486.9 


205.3 856.5 

204.5 866.6 


211.7 

209.3 


977.0 470.6 

993.4 464.5 

1006.4 476.6 

1027.3 480.2 

1043.1 477.7 


2229 3.180 

2269 3.200 

2311 3.222 

2354 3.244 

2399 3.268 


196.3 
257.9 
262.0 

352.4 
38.4 


Table 4 

Computational performance in robust model and robust-MILFP model. 


Model and algorithm Total cost (MM $) Unit cost ($/GGE) a Time(s) 

2354 3.244 109 

2413 3.072 2893 

2414 3.073 (2.984) 36,000 c 

2413 3.072 3033 

2580 3.269 (2.737) 36,000 c 

~(2.786) 36,000 c 

a With the lower bound in the parentheses. 
b Reformulation means reformulation-linearization method. 
c Computational time limit of 10 h is applied to all the algorithms. 


Robust 

Robust-MILFP 

model 


Parametric 

Reformulation 11 

DICOPT 

SBB 

BARON 12.7 


Gap No. of binary 


0% 422 

0% 422 

2.85% 560 
0% 422 

16.27% 422 
422 


No. of continuous 


19.726 
19,729 
20,149 

19.727 
19,727 
19,727 


No. of constraints 


33.835 
33,838 
41,607 

33.836 
33,836 
33,836 


In this robust model, as shown in Fig. 9, 11 preconversion facili¬ 
ties are built, with the capacities ranging from 352 to 648 kton per 
year. The ones in Cook, LaSalle, and Whiteside Counties use the 
technology of fluidized bed reactor, while others use rotating cone 
reactor pyrolysis technology. Eight biorefineries are constructed 
with the capacity ranging from 16 to 93 MM GGE per year, all 
of which use fast pyrolysis followed by hydroprocessing. All the 
four petroleum refineries are retrofitted (Table 2). Note that due to 
the conservatism manner in robust optimization, it tends to have 
tighter demand bounds and lower biomass availability compared 
with the deterministic model. Therefore, it has higher annualized 
cost and the unit cost. In the robust model, more preconversion 
facilities and biofuel refineries are built, and it has higher petroleum 
refinery utilization. 

6.4. Integrated optimization framework 

The last model considered is the robust-MILFP model, which 
integrates the robust counterpart optimization with unit cost 
objective MIFLP model. It consists of the objective function (1), 
the robust counterpart (26)—(29), (31)-(32), and other constraints 
(S2)-(S51), (S54)-(S80) in the supporting information. It can be 
solved by either the tailored algorithms, or the general-purpose 
MINLP solvers. The solution strategy is presented in Fig. 2. The com¬ 
putational results are shown in Table 4. The parametric approach 
returns an optimal solution of $3,072 per GGE within 2893 CPUs. 
The reformulation-linearization approach gets the suboptimal 
solution at the 10 h time limit. Although there is still 2.85% gap 
remaining with $2,984 lower bound, the result of $3,073 per GGE 
is very close to the optimal one. For general MINLP solvers, DICOPT 
still outperforms other solvers. It takes 3033 CPUs to get the opti¬ 
mal solution, which is slightly longer than the parametric approach. 
Within 10 h time limit, SBB returns the result of $3,269 with 16.27% 
gap (lower bound of $2,737). BARON does not return any feasible 
solutions within 10 h. These show that the parametric approach is 
the best choice to this problem. DICOPT is a good alternative and 
can return the optimal solution for this problem although it could 
not guarantee the global optimality for all the models. 

The optimal design of robust-MILFP model is shown in Fig. 10. 
12 preconversion facilities are constructed with the capacity ran¬ 
ging from 297 to 1052 kton per year. Four preconversion facilities 
located in Lee, Pike. Vermilion, and Wayne Counties use the tech¬ 
nology of fluidized bed reactor, while others use rotating cone 
reactor pyrolysis technology. Eight biorefineries are constructed 






































































138 


K. Tong et al. / Computers and Chemical Engineering 68 (2014) 128-139 



with the capacity ranging from 14 to 89 MM GGE per year, all 
of which use fast pyrolysis followed by hydroprocessing. All the 
four petroleum refineries are retrofitted with their full unused 
capacity to upgrade biofuels. Among the four models presented in 
this paper (deterministic, MILFP, robust, and robust-MILFP mod¬ 
els shown in Fig. 6 and Table 2), only in the robust-MILFP models, 
all the remaining capacity in four petroleum refineries are fully uti¬ 
lized. This utilization percentage is much higher than those without 
unit cost objective (the deterministic model and the robust model). 
From the cost breakdown information in Fig. 5, the harvesting cost 
and production cost account for 34% and 30% of the total cost, and 
are the most expensive components. The 0 & M account for the least 
percentage of 7% of the total cost. Although, the models using dif¬ 
ferent objectives could result in different solutions, the cost break¬ 
down will not change too much. The market share information 
is shown in Fig. 7. The petroleum refineries produce 68.9% of the 
biofuels, which is the largest proportion among all the four models. 

7. Conclusions 

This paper addresses the robust design and planning of the 
advanced hydrocarbon biofuel supply chain integrating with 


existing petroleum infrastructure considering the unit cost objec¬ 
tive. Two tailored algorithms are adopted in solving the resulting 
MILFP problems. Demand and supply uncertainties are incor¬ 
porated in the robust optimization framework. Comprehensive 
comparisons between overall cost minimization objective and 
unit cost objective, and between deterministic model and robust 
optimization model are presented respectively. The integrated 
model simultaneously integrates MILFP and robust optimization 
are discussed. Results show the unit cost objective prefers the 
preconversion to petroleum-upgrading pathway, which leads to 
unit cost saving to some extent. Moreover, by controlling the bud¬ 
get parameter r, the robust counterpart programming provides 
a mechanism to control the tradeoff between conservatism and 
economic performance. 
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